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Static potentials for quarkonia at finite temperatures 



We review non-perturbative static potentials commonly used in potential models for quarkonia 
at finite T. Potentials derived from Polyakov loop correlators are shown to be inappropriate for 
this purpose. The qq free energy is physical but has the wrong spatial decay and perturbative 
limit. The so-called singlet free energy is gauge dependent and unphysical. An appropriate 
static real time potential can be defined through a generalisation of pNRQCD to finite T. In 
perturbation theory, its real part reproduces the Debye-screened potential, its imaginary part 
accounts for Landau damping. Possibilities for its non-perturbative evaluation are discussed. 
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1. Introduction 

The properties of quarkonia are believed to provide a useful probe of the QCD plasma 
at high temperatures, in particular for the quark-hadron transition. This expectation was 
originally based on a potential model pQ, in which the linearly confining potential for 
zero temperature gets replaced by a Debye-screened potential at high T. 

Potential models have a long history for the description of quarkonia at zero tempera- 
ture. The basic idea is that for heavy quarks of mass M , which move non-relativistically, 
the binding energy (E — 2M) is small compared to M and can be obtained by solving a 
static Schrodinger equation 



V(r) is the (radially symmetric) potential between the static quark anti-quark pair sep- 
arated by a distance r. Initially V(r) was modelled by the Cornell potential (Coulomb 
plus linear), more recently non-perturbative lattice data are used as input. The crucial 
observation is that the Schrodinger equation follows from an effective theory approach. 
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Starting from QCD, one can use of the scale separation between the heavy quark mass 
M and the binding energy E — 2M, to obtain an effective theory, pNRQCD [2], for 
the low energy dynamics in the confining potential. In this framework, the static po- 
tential appears as a perturbative matching coefficient of the effective theory. Hence, the 
Schrodinger equation can be improved systematically by computing higher order terms 
in the scale hierarchy. Note, that a very successful spectroscopy with ~ 1% accuracy is 
obtained in this way. 

It is tempting to employ this approach also at finite T. Matsui and Satz heuristically 
used the same equation, but with a Debye-screened potential from perturbation theory, 

a 2 Cp c - mo ^ r 

V(r,T) a . (2) 

47r r 

However, there are a number of problems. Firstly, it is not clear if the bound state 
Schrodinger equation can be translated to a finite T many body situation, in a way 
that temperature effects show up only in the potential. Secondly, at finite T there exists 
a variety of non-perturbative potentials, and it is not clear which one constitutes the 
non-perturbative generalisation of Eq. ^ . 



2. Static potentials from the lattice at zero and finite T 



At T = 0, the static potential can be defined non-perturbatively on a euclidean L 3 x N T 
space time lattice. Consider a meson correlation function with an interpolating operator 
^(x)t/(x, y)^(y), where U denotes a straight line gauge string between the quarks. In 
the limit M — ► oo the heavy quarks can be integrated out, taking the correlator to the 
euclidean Wilson loop, 

(VS(x,r)C/(x,y;r)V(y,r)^(x,0)C/(x,y;0)^(y,0)) — > e^ 2MT M/ £ (|x - y|) . (3) 

Inserting a complete set of eigenstates of the Kogut-Susskind Hamiltonian (in temporal 
gauge), the Wilson loop evaluates to (r = |x — y|, U r = U(x, y; 0)) 

We(t,t) = i^|(n|C/ r | TO )|2 e ^„^ e -(^M^„)r (4) 

in 

with E m (r) eigenvalues in the sector with sources, and E n in the sector without. On 
the lattice, T = l/(aN T ), hence T — implies N T — > oo in the second line. Taking 
furthermore the limit r — > oo, the sum is dominated by the lowest energy state. The static 
potential is defined to be the lowest energy of the static quark anti-quark configuration 
at a given separation, V(r) = Ei(r) — Eq- Note that the matrix element with the string 
operator is of no interest here. 

The generalisation to finite T is difficult to interpret because of the finite and short 
temporal extent, N T — \/(aT). Thus, we have to deal with the full superposition Eq. (01, 
to which now also the matrix elements contribute, and the result still depends on r. 

A different definition of the static potential which does generalise to finite T is based 
on the Polyakov loop i(x) = [] T=1 N Uq(x.,t), i.e. a static quark sitting at x and 
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Fig. 1. Static quark anti-quark free energy/potential, Eq. J6jl, for T < T c (left) and T > T c (middle) [5]. 
Right: Free energies for the three channels Eq. J!)}- The solid line is the zero temperature potential [§]• 

propagating in euclidean time through the periodic boundary. It transforms in the adjoint, 
so its trace is gauge invariant. By spectral analysis one establishes that the Polyakov loop 
correlator represents the free energy of a static quark anti-quark pair separated by r [3] , 

e-^ryr = _L (Tri t (x) Tr L(y)) = J_ £ e -*.(r)/T (6) 

c c n 

The energy levels entering this Boltzmann sum are identically the same as the E m (r) 
from the Wilson loop, Eq. Hence, for T — * we recover V(r), cf. [I]. The free energy 
is thus often called a T-dependent potential, V(r,T) = Fq q (r,T). The Polyakov loop 
correlator is readily simulated, with results as in Fig. [T] It gives a linear potential in the 
confined phase, whose string tension reduces with temperature, while in the deconfined 
phase the potential is screened, Fig. [1] Unfortunately, this is not the Debye-screened 
potential we want, as becomes apparent when considering its spatial decay at high T. 
Fitting to 



Fqq _ C ( T ) -m(T)r 

T ' (rT) d 



(J) 



gives d w 1.5 and m = M n ++, i.e. the screening mass corresponds to the lightest, gauge- 

invariant glueball channel [B] . This can already be seen in perturbation theory, where the 
leading term is by two-gluon exchange and thus m = 2mu [Ti- 
lt was thus suggested to decompose the Polyakov loop correlator into channels with 
relative colour singlet and octet orientations of the quark anti-quark pair [3] , 



F q9 (r,T)/T _ 1 -Fi(r,T)/T , § Fs (r,T)/T 

9 9 
W)/T = b TrL i {x)L{y)h 



e -F g9 (r,T)/T _ - e -Fi<T,T)/T + » £ -F a (r,T)/T ^ ^ 
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e- F8(r ' T)/T = J<TrLt(x)TrL(y))-l(TrLt(x)L(y)}. (9) 

Note that the correlators in the singlet and octet channels are gauge dependent, and the 
colour decomposition only holds perturbatively in a fixed gauge. However, in perturba- 
tion theory the singlet channel indeed displays the expected Debye-screened behaviour, 
Fi(T,r) ~ e -m£, ( T ' r /47rr. This has motivated lattice simulations of these correlators 
in fixed Coulomb gauge, with results as in Fig. [T] (right). The three different channels 
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show different r-dependence, and hence lead to different binding energies when used in 
Schrodinger equations. There is a vast literature employing F\ or the corresponding in- 
ternal energy U\ = F± + TSi, and from the solutions trying to reconstruct lattice meson 
correlation functions to check which fits better [9]. 

However, both options are unphysical at a non-perturbative level. To understand this, 
let us start from something physical and consider a meson operator in an octet state, 
O a = i/>(x)[/(x, Xo)T a t/(xo, y)ip(y), with xq the meson's center of mass. In the plasma 
the colour charge can always be neutralised by a gluon. In the correlators for the singlet 
and octet operators, we integrate out the heavy quarks, replacing them by Wilson lines, 



(0(x, y; 0)0* (x, y; N T )) cx (Tr L t (x)C/(x, y; 0)L(y)(7 t (x, y; N T )), 

^J_(TrLt(x)TrL(y)) (10) 

■ (Tr Lt (x) U (x, y ; 0)L (y ) & (x, y ; N T ) } 



(O a (x,y;0)O a t(x,y;iV T ))« 



N C (N? - 1 



We have now arrived at gauge invariant expressions, because we used a gauge string 
between the sources. The singlet correlator corresponds to a periodic Wilson loop which 
wraps around the boundary. The connection to the gauge fixed correlators is readily 
established, replacing the gauge string by gauge fixing functions, J7(x,y) = g~ 1 (x)g(y). 
Thus, in axial gauge t/(x, y) = 1 (and only there) the gauge fixed correlators are identical 
to the gauge invariant ones. 

Next, let us perform the spectral analysis. While indeed the energy eigenvalues in the 
spectral sum are independent of the operators [10] , the full correlators take the form [11] 



c- Fl( - T)/T = ^E^I^(x,y)^(x,y)|n, Q )e-^M/^, 

c n 

e- Fs( - T)/T = ^E^'7l^(x,y)^(x,y)|n, Q )e-^M/ r . (11) 

c n 

The energy levels in the exponents are identically the same in Eqs. (|6I11[) and correspond 
to the familiar gauge invariant static potential at zero temperature and its excitations. 
However, while Eq. ^ is purely a sum of exponentials and thus a true free energy, the 
singlet and octet correlators contain matrix elements which do depend on the operators 
used, thus giving a path/gauge dependent weight to the exponentials contributing to 
Fi,Fg. This is illustrated numerically in Fig. [2] in the low temperature limit, where the 
ground state potential dominates and one can cleanly separate the exponential and the 
matrix elements. The r-dependent structure is entirely in the matrix elements, which 
depend on operators and/or the gauge. 

I do not see how this is evaded by applying smearing techniques, as recently suggested 
in [12) . These authors replace the spatial string swith a smeared object in order to 
increase the overlap with the ground state, i.e. to get the ground state matrix element 
close to one. However, most smearing changes the expectation values of correlators, thus 
destroying their mutual relations, Eq. ([9]). Secondly, smearing increases the weight of 
the lower energy states at the cost of the higher ones, and thus undoes the effect of 
finite temperature in a procedure dependent way. Finally, if one could get all matrix 
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Fig. 2. Left: Polyakov loop correlators, Eq. JS}, for r/a = 1 in the case of 3d SU(2) in the low temperature 
limit. All decay with the same ground state exponential. Right: The corresponding matrix elements in 
the three channels introduces operator dependent r-dependence, except for the average channel 

elements equal to one, the different channels would simply be equal, up to the trivial 
colour coefficients, with no additional information. 

To summarise, since the spectral information contained in the average and gauge fixed 
singlet and octet channels is the same, we must conclude that any difference between 
those correlators is entirely gauge dependent and thus unphysical, and so are all binding 
energies calculated in F\ oxU\. 



3. A real time static potential for finite T quarkonia 



Progress was made recently by generalising the effective theory approach quarkonium 
physics at T = 0, namely pNRQCD, to finite temperatures [13114115] . as reviewed at this 
conference [16]. The analysis is performed in a perturbativc setting in Minkowski time. 
Just as at zero temperature, the static potential then appears as a matching coefficient 
in the effective theory after the heavy modes have been integrated out. The relevant 
correlation function is the quarkonium correlator in real time, but evaluated as a thermal 
expectation value. Not surprisingly, after integrating out the static quarks, the correlator 
is proportional to a Wilson loop in Minkowski time, We(H, r). Of course, the expectation 
value implied in We is now a thermal one, i.e. N T is finite for fixed lattice spacing. Hence 
we need the analytic continuation of the double spectral sum in Eq. {3J . From the effective 
theory it is easy to see that this correlator obeys a real time evolution equation 

[id t -V > (t ) r)]W E (it,r)=0. (12) 

This represents the desired Schrodinger equation for quarkonia in the plasma, and defines 
the relevant real time dependent potential. The required scale hierarchy for this equation 
to be valid is g 2 M < T < gM. Furthermore, for non-relativistic bound states p <C E, 
hence we need Or, i.e. the static pontential is obtained in the long time limit V(oo, r). 

Eq. (|12p may be also be viewed as a non-perturbative definition of the potential of in- 
terest via a correlation function, just as was the case for the zero temperature potential. 
Unfortunately, this one is defined for Minkowski time and thus requires analytic contin- 
uation, i.e. it cannot be evaluated directly from euclidean lattice simulations. However, 
a first impression about this object can be gained form HTL-resummed perturbation 
theory, for which the leading order result is 
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V>(oo,r) 



9 2 C F 



4tt 

with 4>{x) = 2 



dz z 



exp(-m D r) 



sin(za;) 



47T 



(z 2 + l) 2 



(m D r) , 



(13) 



The most striking feature of this potential is that it is complex, contrary to the free 
energies discussed before. The real part features the expected Debye-screened potential. 
The imaginary part is due to Landau damping and must necessarily be there for a correct 
effective description of the plasma dynamics. Its derivation and properties are discussed 
in more detail in [13114116] , 

What are the corrections to this potential? Firstly, there are corrections from HTL- 
resummed perturbation theory of the order g 2 T/A, where A is the UV cut-off, with a 
calculable coefficient. Here, we are interested in the non-perturbative corrections from 
infrared modes ~ g 2 T/m mag . These are due to the soft colour magnetic modes m mag ~ 
g 2 T, and thus cannot be calculated in perturbation theory. 

However, one can calculate these non-perturbative corrections by classical lattice sim- 
ulations in that sector of the theory, which has high occupation numbers and is well 
represented by a classical approximation. To identify this sector it is instructive to take 
the limit H — + in our perturbative result Eq. (fTS"]) first. To this end, H needs to be 
reinstated by the replacements g 2 — > g 2 h, 1/T — » H/T, leading to 



lim V>(oo, r) 



ig 2 TC F 

Air 



cf>(m D r). 



(14) 



Thus, only the imaginary part survives in the continuum limit. This is easy to understand 
since the long range physics of Landau damping is dominated by classical fields, e.g. in 
scalar field dynamics, whereas the binding is a generic quantum effect, cf. the hydrogen 
problem. Thus, we can evaluate the non-perturbative infrared effects for the imaginary 
part of the potential. 



4. Imaginary part from classical lattice simulations 



This has been done in [17], following the technical setup that was also used for the 
evaluation of the sphaleron rate in the electroweak theory [TH]. In order to perform 
real time simulations one reformulates the theory in a Hamiltonian approach. Fixing 
temporal gauge Uq — 0, the conjugate field operators are the links and the electric fields 
defined by U(x,t) — iEi(x,t)Ui(x,t). Full gauge invariance is restored by imposing the 

Gauss constraint G(x) = J^i Ei( x ) ~ U-i(x)Ei(x — zjf/l^i) — j°(x) — 0. A thermal 
distribution at some initial time is generated by the partition function 



Z=[vU i VE, S(G)e-? H , H = ±- £ Re Tr(l - U {j ) + \ Tr(£ 2 ) 

x i<j 



, (15) 



where Uij is the plaquette. The distribution is then evolved in real time by the classical 
equations of motion, 
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Fig. 3. The imaginary part of the real-time static potential from the classical simulation (left) and from 
resummed perturbation theory (right) 17 . 



U i (x)=iE i (x)U i (x), Ei 



E?T a , Ef(x) = -2ImTr[T Q - ]T U l3 {x)] . (16) 



This procedure can be improved by taking into account quantum corrections by using a 
HTL-rcsummcd effective theory. It will generate a source term due to the hard particles 
in the plasma, which modifies the Hamiltonian as well as the equations of motion by 
coupling the classical fields to the quantum effects of the hard particles. That approach 
was used in [T5] and is easily adapted to the present problem [T7] . 

One can now calculate the classical thermal expectation value of the real time Wilson 
loop, which in temporal gauge reduces to a correlation function of the spatial string, 



W E ,cl(it,r) = (f/ t (x,y;t)[/(x,y;0)). 



(17) 



The time dependent potential extracted via Eq. (JT^J) is shown in Fig. [3J where it is 
compared with the result obtained from resummed classical lattice perturbation theory. 
We observe complete qualitative agreement. Note in particular that a finite imaginary 
part survives also non-perturbatively in the long time limit, when fluctuations have died 
out. Comparing the values one finds that the non-perturbative effects make the imaginary 
part more negative, i.e. increase Landau damping. 



5. Non-perturbative real time potential? 



It is clearly desirable to go beyond the classical approximation and construct an opera- 
tor from which the whole quantum potential, including its real part, can be extracted. Of 
course, a full quantum computation of the real time Wilson loop is impossible, just as the 
calculation of any other real time correlation function. The whole point of the potential 
approach is to bypass the need for such correlators. In our case, what is needed is the 
static potential in the infinite time limit, which is clearly less information than having 
to know the full time dependence. In terms of correlators, the information we need is 

lim W E (it,r), lim d t W E (it,r). (18) 

f— »CO t — >oo 
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At least in principle, these limits ought tho be representable by Euclidean operators, the 
challenge is to construct those in practice. 

6. Conclusions 

We have argued that many potential models used for the description of quarkonia 
at finite T have significant flaws. The connection between the Schrodinger equation to 
the underlying quantum field theory is unclear, and lattice potentials extracted from 
Polyakov loop correlators, which are typically used as input for those models, are the 
wrong quantities for this purpose. The average free energy is gauge invariant and well 
defined, but in its perturbative limit does not reduce to the Debye-screened potential. 
The so-called singlet potential is gauge dependent and therefore unphysical. 

These problems can be overcome by using an effective field theory obtained by inte- 
grating out the heavy quarks, which is pNRQCD generalised to finite temperatures. The 
resulting Schrodinger equation is the real time evolution equation for a quarkonium cor- 
relator, and the static potential in this equation is a matching coefficient in the effective 
theory. To leading order in HTL-resummed perturbation theory, this potential is complex, 
its real part showing the correct Debye-screened behaviour and its imaginary part reflect- 
ing Landau damping. In the classical limit, only the imaginary part survives. This part 
can be calculated non-perturbatively with classical lattice simulations in real time. The 
result agrees in all qualitative features with the HTL result with slightly strengthened 
damping. It is now important to search for a lattice operator that represents the real part. 
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